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Abstract 

In this work, we calculate the /i c ( 1 Pi) production rate at the LHC to leading order of the strong 
coupling constant, for both color-singlet and -octet mechanisms. Numerical results show that a 
considerable number of h c events with moderate transverse momentum px will be produced in the 
early run of the LHC, which will supply a good opportunity to further study the nature of this 
P-wave spin-singlet charmonium state. 

PACS number(s): 13.85.Ni, 14.40.Lb, 12.39.Hg, 12.38.Bx. 
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Since the first charmonium, the J/ip, was discovered thirty years ago, much effort has 
been made to explore it and its higher excited states with both theory and experiment. These 
studies have provided deep insights into the heavy quark-antiquark strong interaction, or, in 
other words, the application of quantum chromodynamics(QCD). Although much progress 
has been made, there are still many unsolved problems left in the study of quarkonium 
physics. For instance, in the charmonium sector, the cc mass spectrum of the naive quark 
model prediction has not been completely confirmed experimentally yet. Below the open 
charm threshold, all expected charmonia have been identified in recent years, but experi- 
mental measurements of the physical natures of rf c and /i c ( 1 Pi) are quite limited. The spin 
singlet states of heavy quarkonia pose an experimental challenge because they are not pop- 
ulated at lepton colliders. In hadron-hadron collision, the X P\ state can be formed directly 
in many ways. The goal of this work is to analyze the possibility of detecting /i c ( 1 P 1 ) at the 
LHC. 

h c is the ground state of the P-wave spin-singlet in the charmonium family. According 
to the QCD-based potential model prediction, to leading order of the spin-spin interaction 
the hyperfine splitting AM h/ (M( 1 P 1 ) - M( 3 Pj)) should be zero. Here, the spin- weighted 
average mass of P-wave triplet states M( 3 Pj) = (M Xo ( 3 P ) + 3M Xo ( 3 P!) + 5M X0 ( 3 P 2 ))/9 = 
3525.30±0.04MeV and higher order corrections to the hyperfine splitting should be less than 

nnn 

1 MeV pj, 12J, 13[. In 1995, the h c signature, at about 3526 MeV, was first observed in the 
channel of h c — > J/ipir by the E760 Collaboration at the Fermilab 4j. Although this result 
was not confirmed by E835, which succeeded E760 with significantly higher statistics, the 
E835 Collaboration reported that they observed evidence of h c via the h c — > rjc'-y process and 



obtained a resonance mass of 3525.8 ± 0.2 ± 0.2 MeV 



5]. In the electron-positron collision, 



the CLEO Collaboration reported that they measured the mass of h c at 3525.28±0.19±0.12 



MeV via the decay of ip(2S) -> n°h c followed by h c -> r/ c7 at CESR Q, H, E], while the Belle 
Collaboration did not observe significant signal in the decay of B — > h c K ± 9J. For a more 
detailed theoretical description of h c and experimental progress in this respect, readers are 



referred to reviews 



10 



11 



12|. 



To obtain more knowledge of the nature of h c , a key point for experimentalists is to 
obtain enough h c event data. The Large Hadron Collider (LHC) will be operational this 
year, which may supply a good opportunity to study quarkonium physics, including h c . 



2 



fa) (b) (c) 

FIG. 1: Typical Feynman diagrams of h c production in the extrinsic charm- induced process g+c — > 
/i c ( 1 p| 1 ^) + c in the color-singlet scheme. 

With a luminosity of about 10 32 ~ 10 34 cm~ 2 s™ 1 and a center of mass energy of 10 ~ 14 
TeV, the LHC will produce copious charmonium data, which in principle will enable people 
to measure the h c state more precisely. In the following we evaluate the h c production rate 
at the LHC . 



It is well-known that historically the so-called color-singlet model(CSM) [13l. |14|. Il5l. |16, 



17l | played a major role in the study of quarkonium physics and had great success in many 
respects. However, it failed to explain the Fermilab Tevatron data of charmonium large 
transverse momentum production. Hence, the color-octet mechanism (COM) was proposed 



and employed 



18 



191 ]. which is based on a solid framework, the nonrelativity QCD(NRQCD) 
20j . The effective theory of NRQCD is widely accepted nowadays, although the validity of 
applying it to the charmonium phenomeno logical study is in some sense still vague. In the 
following calculation, nevertheless, both color-singlet and -octet contributions will be taken 
into account. 

The differential cross section for h c hadroproduction is formulated in a standard way, 
•j-^(PP -> h c + X) = J dx a dyf a/p (x a )f b/p (x b ) ^~Z~^^ ^( a + b -> hc + X ) > (*) 

a,b 

where f a / p and f b / p denote the parton densities; s, t, and u are Mandelstam variables at the 
parton level; y stands for the rapidity of produced h c ; xt = with mx = \/ M 2 + p^; 
and the capital v^S* and M denote the total energy of incident beam and the mass of h c , 
respectively. 

To leading order and with moderate transverse momentum, the dominant partonic sub- 
processes for h c hadroproduction evidently include 

g + g _> h c CS [ * ] ) + g , (2) 
9 + q(q) - h c CS [ * ] ) + q(q) , (3) 



q + q 
9 + 9 
9 + c(c) 



+ g 



hCPl 1 ') + 9 , 
h c CPl 1] ) + c(c) , 



(4) 
(5) 
(6) 



where the first three represent the h c production processes in the color-octet scheme, while 
the last two are through CSM. The process (jSj) is an "extrinsic charm" one, and its im- 



portance in charmonium hadroproduction was exhibited in Refs. 



21 



22 



23j |. To the lowest 



order of the strong coupling constant, expressions for the partonic differential cross section 



da/dt of processes (J2]) to (j5J) were obtained in several previous studies |24j, |25|, |26|, |27|, |28| . 
whereas the analytic expression for process g + c — > /i c ( 1 P 1 ' 1 ^) + c, as schematically shown in 
Figure 1, is still absent in the literature. It is worth mentioning that the "extrinsic charm" 
induced process c + c — > h c + g is omitted in our calculation since its numerical contribution 
is negligibly small. 

For process ([6]), we commence with the calculation of the partonic process g + c — > 
(cc) + c, then project the cc matrix element onto the color-singlet 1 P] 1 ' state. In calculating 
processes involving P-wave heavy quarkonium to leading order accuracy in relativistic expan- 
sion, one must expand the amplitude to the second order in powers of the relative momentum 
between the constituents of heavy quarkonium since the first order term gives no contribu- 
tion. After taking the non-relativistic limit, it is then legitimate to take p c = p 5 = P/2, 
one half of the charmonium momentum produced. For the outgoing h c , one can employ the 
following projection operator: 



v(Pe) u(p c ) 



1 
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P 

~ m c) 75 (y- 4+ m c 



(7) 



2V^ m c 

where q is the relative momentum between two charm quarks, N c = 3, and l c represents 
the unit color matrix. By writing the projector in a matter similar to (CO), it is understood 
that M = 2m c has been implicitly assumed. 

After following the procedures mentioned above, it is straightforward to calculate the 
process (EJ), and the analytic result reads 

12 / 9t 



da l6a 3 s n\R'( 
dt 27m c (s — m 2 ) 2 



's — ml) 2 ml 
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FIG. 2: The h c production rates as a function of the transverse momentum lower bound pt in pp 
collision at the center-of-mass energy y/S=U TeV. The left dia gram demonstrates the integrated 
cross-sections of h c production via processes ([2]) to (|6|) shown as lines a to e, respectively. The 
solid line in the right diagram represents the yield from the color-octet scheme, and the dashed 
line represents the yield from the color-singlet scheme. 
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Here, the nonperturbative parameter, i?^ (0), is the derivative of the Schrodinger radial wave 
function at the origin for h c , which can be either inferred from phenomenological potential 
models or extracted from experimental data. 

In our numerical evaluation, the input parameters are taken as follows: v^? = 14TeV, 
m c = M/2 = 1.78GeV, the value of the color-singlet matrix element (0|C?i c ( 1 Pi)|0) = 0.32 



GeV 5 



29j, the value of the color-octet matrix element (OlOg^ 1 ^)^) = 9.8 x KT 3 GeV 3 
28l |. and the pseudorapidity cut \rj(h c )\ < 2.2 is enforced according to the LHC experimental 
environment. In the calculation, the typical energy scale is set to be at mx = ^/M 2 + p T ; the 
strong coupling constant a s is running with transverse momentum. Both renormalization 
and factorization scales are evolved to the same point mx, and the CTEQ5L |30|] parton 
distribution function is employed. In Eq.pl), the relation between the NRQCD matrix 
element and the derivative of the Schrodinger radial wave function at the origin for the X P\ 
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TABLE I: h c production rates with various transverse momentum lower bounds at the center-of- 
mass energy = 14TeV and an integrated luminosity of 10 fb _1 are presented. Taking into 
account the three main decay chains of h c , i.e. 1) h c — ► 7ToJ/V' — ► [J- + 7 i_ 77> 2) h c — > rjc'j — > p|ry 
and 3) /i c — ► 7^7 — ► 777, the final experimentally detectable event numbers are given. 





color-singlct event 


color-singlct event without charm sea effect 


color-octet event 


VTcut 


5 GeV 


10 GeV 


20 GeV 


30 GeV 


5 GeV 


10 GeV 


20 GeV 


30 GeV 


5 GeV 


10 GeV 


20 GeV 


30 GeV 


total || 1.65 X 10 8 | 4.32 x 10 6 | 8.14 x 10 4 | 7.57 x 10 3 | 8.41 x 10 7 | 1.41 x 10 6 | 1.02 x 10 4 |4.70 X 10 2 | 3.78 X 10 9 1 1.56 x 10 s | 3.67 x 10 6 | 8.54 x 10 5 


Chains 


4.94 X 10 4 


1.30 x 10 3 


2.44 x 10 


2.27 


2.52 X 10 4 


4.22 x 10 2 


3.06 


0.14 


1.13 X 10 6 


4.68 x 10 4 


1.10 x 10 3 


1.06 X 10 2 


Chains 


1.07 x 10 s 


2.81 X 10 3 


5.29 x 10 


4.92 


5.47 x 10 4 


9.14 x 10 2 


6.64 


0.31 


2.45 X 10 6 


1.01 X 10 5 


2.38 X 10 3 


2.30 X 10 2 


Chains 


1.97 x 10 4 


5.19 x 10 2 


9.76 


0.91 


1.01 X 10 4 


1.69 X 10 2 


1.23 


0.06 


4.53 X 10 5 


1.87 x 10 4 


4.40 X 10 2 


4.24 x 10 



state, i.e., 

I#(0)| = ^(01^(^)10) , (9) 

is adopted. Note that among the inputs, the charm quark mass m c is taken to be one half 
of the h c mass for simplicity, i.e., the constituent quark mass, which we find may increase 
the final result by some 30% from that found when taking m c to be 1.5 GeV. 

The numerical results of the integrated cross section for different pt lower bounds are 
given in Figure [2j From the figure, it can be found that the contribution from COM is about 
two orders of magnitude larger than that from CSM in almost every transverse momentum 
region. Among the three color-octet processes, the contribution from process ([2]) dominates 
over the other two. Of the two color-singlet processes, the yield from process (EI) overshoots 
that from process ([5]) in the large transverse momentum region, in spite of the suppression of 
the extrinsic charm distribution. Because of the big gap between the yields from the color- 
singlet and color-octet, one result of this calculation is that the experimental measurement 
may tell whether the color-octet estimate of h c production is reliable or not. 

In experiment the h c can be reconstructed from its three dominant decay modes, which 

are 

h c -> 7T J/tp -> yuV~77 , (10) 

h c -> VcJ -> PP1 , (11) 
h c -> Vd -> 777 • (12) 

Of these decay chains, J/ip decays into with a branching ratio of 6% [3jJ, 7r° almost 
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completely decays into 77, and rj f . decays into pp with a branching fraction of 0.13% and 
into 77 with a ratio of 0.024% [31]. The branch fractions of h r — » J film and h c — > rj c ^ are 



theoretically estimated to be about 0.5% [32J and 50% |33j, ]34|, [35j, ]36|] , respectively. For the 
h c —> J/ip7T° process, although the 7r°s produced are energetic, their decays to two photons 
can be well resolved when the ir° momentum is less than 40 GeV [33] . Considering the 
decay rates of h c to these experimentally measurable modes, in Tabled] we present the event 
numbers of the decay chains (II 01) - (II 21) with different transverse momentum lower bounds and 
in the LHC experiment environment, that is, a 14 TeV colliding energy, a 10 fb _1 integrated 
luminosity and a pseudo-rapidity cut \r](h c )\ < 2.2. From the table we see that even the 
h c produced with a lower transverse momentum bound of 10 GeV, in which region the 
experimental detection efficiency becomes high, there will be millions of events coming out 
in its three dominant decay modes, from both the color-singlet and -octet schemes. In the 
table, we also present the color-singlet contribution without the charm sea effects. One may 
find that the charm sea-induced process contributes at least half of the total color-singlet 
yield with various transverse momentum lower bounds. 

In conclusion, we have evaluated the h c direct production rate at the LHC, where the 
h c indirect yields are much less than the direct ones according to a similar analysis for h c 
production at HERA-b 271 ] . Our calculation is performed to leading order of the strong 
coupling constant a s and to second order in the relative velocity v 2 expansion. Both color- 
singlet and -octet production schemes are taken into account in this work. We find that there 
will be enough h c yields at the LHC for a precise measurement on the nature of this P-wave 
spin singlet. Although as usual the high order corrections may induce some uncertainties 
in the calculation, as an order-of-magnitude estimate our results should hold. Due to the 
large discrepancy between predictions from the color-singlet and color-octet schemes, the 
experimental measurement of the h c production rate at the LHC may tell to what degree 
the color-octet mechanism plays a role in charmonium production as well. 



Finally, as we were studying this issue, there appeared a similar work on the web 38]. 



The main difference between this work and Ref. 



contribution process (J6]). Since the necessary definitions in several places of Ref. 



is the inclusion of the extrinsic charm 

are not 



clear, it is hard to make a direct comparison of our results with those given in the reference. 
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